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Abstract. 

High resolution, multi- wavelength maps of a sizeable set of nearby galaxies have made it 
possible to study how the surface densities of Hi, H2 and star formation rate (Ehi,Eh2,Esfr) 
relate on scales of a few hundred parsecs. At these scales, individual galaxy disks are comfortably 
resolved, making it possible to assess gas-SFR relations with respect to environment within 
galaxies. Eh2, traced by CO intensity, shows a strong correlation with Esfr and the ratio 
between these two quantities, the molecular gas depletion time, appears to be constant at about 
2 Gyr in large spiral galaxies. Within the star-forming disks of galaxies, Esfr shows almost no 
correlation with Em. In the outer parts of galaxies, however, Esfr does scale with Em, though 
with large scatter. Combining data from these different environments yields a distribution with 
multiple regimes in S gas — Esfr space. If the underlying assumptions to convert observables 
to physical quantities are matched, even combined datasets based on different SFR tracers, 
methodologies and spatial scales occupy a well define locus in E gas — Esfr space. 
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1. Introduction and the Global Star Formation Law 

Great progress has been made towards an understanding of star formation (SF) on 
small scales in the Milky Way, but many open questions remain about its connection 
to large scale processes: what sets where SF occurs in galaxies and how efficiently gas 
is converted into stars? How important are global, galaxy-scale environmental parame- 
ters as opposed to small-scale properties of the interstellar medium (ISM)? What is the 
role of feedback in regulating SF? To address such questions, theoretical modeling and 
simulations need to be constrained by comprehensive observations. 
t-H . Both observations and theory have focused on the relationship between the star forma- 
tion rate (SFR) and the gas density, for which a tight power-law relationship was observed 



• ■ *h . in a large number of galaxies by Kennicutt(1998) Such a relationship was first suggested 
many decades ago by Schmidt(1959)| who studied the distributions of atomic gas and 



stars in the Galaxy. Over the following decades, similar studies targeted individual Local 



Group galaxies, e.g., M33 (|Madore et al.(1974)| |Newton(1980)| , the Large Magellani c 
Cloud (|Tosa fc Hamajima(1975)|), and the Small Magellanic Cloud (|Sanduleak(1969)[). 



Kennicutt(1989) carried out the first comprehensive extragalactic study targeting a large 



sample of nearby galaxies and Kennicutt(1998) followed up this work, focusing on mea- 



surements averaged across galaxy disks. In a sample of 97 nearby normal and starburst 
galaxies, he found a close correlation between the galaxy-average total gas surface den- 
sity (Sgas = Shi + SH2) and the galaxy-average SFR surface density (Esfr)- Following 
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Figure 1. Left: Kennicutt(1998) found a strong correlation over many orders-of-magnitude 
between global averages of Esfr and E gas = Ehi + Em f° r a large sample of nearby normal 
and starburst galaxies. He derived a power law index TV ~ 1.40, implying more efficient SF 
for galaxies with higher average g as columns. Middle a nd Right: Surface densities of gas and 
SFR measured in radial profiles by Wong & Blitz(2002) for two exemplar y nearby spirals. The 
middle panel shows Ssfr versus Ehi, the right panel Esfr versus Eh2- Wong & Blitz(2002) 
found no correlation between atomic gas and SFR, whereas molecular gas and SFR scale with 
one another. 



this work, it has become standard to study the relationship between gas and star forma- 
tion via surface densities, which are observationally more easily accessible than volume 
densities. 



Kennicutt(1998) found Ssfr — Ax S^ s , with intercept A and power law index N — a 
relationship that is variously referred to as the "star formation law," "Schmidt-Kennicutt 
law," or "Schmidt Law." Kennicutt(1998) derived N rj 1.40. Because the ratio SsFR/S gas 
describes how efficiently gas is converted into stars (and is thus often referred to as the 
star formation efficiency, SFE), this super- linear power law index implies that systems 
with higher average gas surface densities more efficiently convert gas into stars (left panel, 
Figure [lj. This measured value is close to N — 1.5, which is expected if the free-fall time 
in a fixed scale height gas disk is the governing timescale for SF on large scales. Other 
studies working with disk-averaged, global measurements found N to be in the range of 
~0.9 



1.7 (e.g., Buat et al.(1989)||Buat(1992)| peharveng et al.(1994)| 



The availability of high-resolution maps of CO emission, the standard tracer of molecu- 



lar gas, made it possible to follow up the work of Kennicutt(1998) with studies focusing on 
azimuthally-averaged radial profiles of gas and SF. Resolving galaxies in this way makes 
it possible to look at how gas and SF relate within individual galaxy disks, opening up a 
wide range of environmental factors to explore. | Wong fc Blitz(2002)|us ed BIMA SONG 
data ( |Helfer et al.(2003)| ) to study 6 nearby spirals, |Boissier et al.(2003)] explored a larger 
sample of nearby spirals, Heyer et al.(2004) studied the Local Group galaxy M33, and 



Schuster et al.(2007) explored the gas-SF relation in M51. These studies derived power 
law indices in the range N rs 1 — 3, leaving it unclear whether a single relation relates 
gas and SF when galaxy disks are spatially resolved. Further disagreement centered on 
the relationship of SF to different types of gas — Hi, H2, and total gas. Intuitively, one 
might expect a stronger correlation between SF and the cold, molecular phase, rather 



than the atomic phase. Wong & Blitz(2002) indeed found a much stronger correlation of 
Ssfr with the molecular gas, Sh2 (compare Figured]). However, Kennicutt(1998) and 



Schuster et al.(2007) both found a better correlation of Ssfr with the total gas, S g 



than with Sjj2- 
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2. Recent Advances: Gas and Star Formation on sub-kiloparsec 
Scales 

One reason that different studies returned such different results were the wide range 
of SFR tracers employed in the various analyses. Furthermore, different studies em- 
ployed widely varying methods to correct the observed UV and Ha intensities for the 
effects of extinction by dust. Because the correction factor is usually > 2, the adopted 
methodology makes a large difference. A large step forward in this field came from the 
Spitzer space telescope. As part of SINGS (Spitzer Infrared Nearby Galaxies Survey, 
Kennicutt et al.(2003)), Spitzer obtained high- resolution IR maps of a large sample of 
nearby galaxies. Calzetti et al.(2005) and Calzetti et al.(2007) demonstrated the utility 
of these maps to trace recently formed stars obscured on small scales, particularly when 
used in combination with Ha — a tracer of unobscured star formation. 



At the same time the VLA large program THINGS (Walter et al.(2008) I obtained the 
first large set of high resolution, high sensitivity 21-cm line maps for the same sample of 
galaxies. Following shortly thereafter, the IRAM large program HERACLES mapped CO 



emission for an overlapping sample of nearby galaxies (first maps in Leroy et al.(2009j| ). 
The result was, for the first time, a matched set of sensitive, high spatial resolution 
maps of atomic gas, molecular gas, embedded and unobscured star formation for a large 
sample of nearby galaxies. The resolution of the maps allowed hundreds of independent 
measurements per galaxy, leading to significantly improved statistics and the ability to 
carefully isolate regions with specific physical conditions. 



Bigiel et al.(2008) combined these data to compare Hi, H 2 , and SFR across a sample 



of 7 nearby spiral galaxies. Figure [5] shows the results of this analysis: the left panel shows 
Shi, the middle panel Sh2, and the right panel S gas = Shi + Sh2 versus Ssfr (derived 
from a combination of far UV and 24/im emission). Galex far UV emission was chosen to 
trace the recent, unobscured SF because of the low background in the FUV channel and 
the large field-of-view of the GALEX satellite. In these plots, H I and H2 show distinct 
behaviors: the atomic gas shows no clear correlation with the SFR, whereas the molecular 
gas exhibits a strong correlation. As a result, the composite total gas-SFR relation is more 
complex than a single power law. In the combined (total gas) plot in the right panel, one 
can clearly distinguish where the ISM is H 1 dominated (low gas columns, steep relation) 
form where it is H2 dominated (high gas columns, roughly linear correlation). 

If a power law is fit to the molecular gas distribution in the middle panel, one obtains 
N « 1.0. This can be restated as a constant ratio Ssfr/Sh2, which means that on average 



each parcel of H2 forms stars at the same rate. Leroy et al.(2008) searched for correlations 



between Ssfr/Sh2 and a number of environmental variables — ISM pressure, dynamical 
time, galactocentric radius, stellar and gas surface density — and found little or no 
variation across the disks of 12 nearby spirals. On the other hand many of these same 
environemntal variables do correlate strongly with the H2-to-Hi ratio. The combined 
conclusion of these two studies was that the average depletion time in the molecular gas 
(i.e., Sh2/Ssfr) of nearby spirals is fairly constant at ~ 2.0 Gyr, while the abundance 
of molecular gas is a strong function of environment. 



Blanc et al.(2009) carried out a similar experiment to Bigiel et al.(2008) They sampled 



the inner part of M51 with 170 pc diameter apertures and estimated local SFR surface 
densities from Ha emission. Their integral field unit observations allowed for accurate 
estimates of internal extinction and corrections for contamination by the AGN and diffuse 
ionized gas. They used a Monte-Carlo approach to incorporate upper limits into their 



power law fit. Their results are in good agreement with Bigiel et al.(2008) regarding M51 



in particular as well as regarding the general conclusions they reached: a virtual absence 
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Figure 2. Esfr versus £ht (left), Eh2 (middle) and E gas (right panel) for pixel- by-pixel data 
from 7 nearby spirals at 750 pc resolution. The contours represent the density of sampling points 
(pixels), where darker colors indicate a higher density. The Hi distribution saturates at about 
10 Mq pc -2 and shows no correlation with the SFR (left panel). Gas in excess of this surface 
density is predominently molecular. The middle panel illustrates the correlation between H2 and 
SFR, which can be described by a power law with slope AT w 1.0. This implies a constant H2 
depletion time of ~ 2 Gyr. The total gas plot (right panel) illustrates the different behavior of 
Hi and H2 dominated ISM at low and high gas column densities, respectively. 



of correlation between SFR and Hi, a strong correlation between SFR and H2 and a 
molecular gas depletion time that shows little variation with molecular gas column. 

Recently, Rahman et al.(2010) explored the impact of different SFR tracers and the 
role of possible contributions from diffuse emission and different sampling and fitting 
strategies on the relationship between Ejj2 and Ssfr- They found that the SFR derived 
for low surface brightness regions is extremely sensitive to the underlying assumptions, 
but that at high surface brightness the result of a roughly constant H2 depletion time is 
robust. 

Even more recently, Schruba et al. (in prep.) combined the HERACLES and THINGS 
data to coherently average CO spectra as a function of radius. With this approach they 
are able to trace molecular gas out to 1.2 r25, allowing them to study the relation between 
H2 and SFR where Hi dominates the ISM. They demonstrate that the tight correlation 
between H2 and SFR crosses seamlessly into the Hi-dominated outer disk (left panel 
Figure El). 

With so much effort expended measuring SFRs and gas densities over the years, it 
is interesting to ask whether the literature largely agrees. The right panel in Figure [3] 



shows a collection of literature measurements along with the data from Bigiel et al. (2008) 



After matching underlying assumptions about how to derive physical quantities from the 
observables, the literature data populate a well-defined locus in £sFR _ E gas space. 

With some consensus emerging on the broad distribution of data in SFR-H2 space, at- 



tention is turning to the origin of the intrinsic scatter in the SFR-H2 ratio. Schruba et al.(2010) 



looked at this as a function of spatial scale in M33 and showed that scatter in the 
CO-to-Hev ratio increases dramatically once a resolution element contains only a single 
star- forming region (i.e., Hn region or giant molecular cloud). This occurs at scales of 
~ 150 pc in M33 but should be a function of the environment studied. They interpreted 
this finding as a result of the evolution of these regions, so that information on the life 
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Figure 3. Left: Esfr versus Eh2 from Schruba et al. (in prep.). They apply a stacking analysis 
to the HERACLES CO data to probe the H2-SFR relation far into the regime where Ehi > Eh2- 
The correlation between Esfr and Eh2 extends smoothly into the H i-dominated parts of galaxies 
out to 1.2 r25. Righ t: Comparison between different datasets using different methodologies from 
Bigiel et al.(2008) The contours are identical to the right panel in Figure [2] and the overplotted 
symbols come from studies using a variety of SFR tracers and methodologies. All datasets have 
been adjusted to match the same set of assumptions when converting observables to physical 
quantities. The composite sample occupies a well-defined locus in Esfr — E gas space. 



cycle of giant molecular clouds is embedded in the Esfr-Eh2 relation, especially at high 
resolutions. 

Another avenue of investigation is the role of the host galaxy in driving the scatter in 



the SFR-H2 ratio. Leroy et al.(2008) and Bigiel et al.(2008) found little environmental 



dependence of this quantity in large spiral galaxies, but with the completion of the 
HERACLES survey (47 galaxies spanning from low-mass dwarfs to large spirals) we can 
now apply a similar analysis to a much wider sample of galaxies. Bigiel et al. (in prep.) 
and Schruba et al. (in prep.) explore how scatter in the SFR-to-H 2 ratio breaks into 
scatter within galaxies and scatter among galaxies. They find that scatter among galaxies 
dominates the relation, with a clear trend evident so that less massive, more metal poor, 



later type galaxies show systematically higher ratios of SFR-to-H2. Young et al.(1996) 
saw this trend using integrated FIR-to-CO ratios. These new analyses show that once 
these galaxy-to-galaxy variations are removed, galaxies exhibit a very tight internal H2- 
SFR relation. 



3. Scaling Relations beyond the Optical Disks 

Hi maps reveal atomic gas out to many optical radii in spiral galaxies and over the 
past few years, GALEX UV observations have revealed widespread SF in the outer parts 
of many galaxies (e.g., |Thilker et al.(2005)||Gil de Paz et al.(2007)||Bigiel et al.(2010b)| ). 
These outer disks have fewer heavy elements, less dust, and lower stellar and gas surface 
densities than the inner parts of spiral galaxies. Contrasting the gas-SFR relationship in 
outer disks with that in the inner parts of normal galaxies gives a chance to assess the 
impact of these parameters on SF. 

In the left panel of Figure 2] we show the results of a pixel- by-pixel analysis of outer 
galaxy disks: the open contours show Esfr versus Egas for a sample of 17 spiral galax- 
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Figure 4. Left: Esfr versus S Eaa for the ou ter (open contours) and inner (filled contours) parts of 



nearby spiral galaxies (Bigiel et al. (2010a) I. The combined distribution reveals multiple regimes: 
at large radii and low surface densities, Esfr scales with E gas ~ Shi. Over most of the area in 
disks, Esfr is a very steep function of E gas , with the kb-to-H I ratio being the key determinant 
of Esfr. At high column densities, the gas is predominantly molecular and correlates well with 
Esfr. At even higher column densities, a steepening of this relation, meaning an increasing 
efficiency of SF, may accompany the transition from galaxy disks to starbursts. Right: Hi and 
far UV radial profiles for M83 out to almost 4 optical radii 125. The inferred Hi depletion time is 
about a Hubble time at large radii, much longer than the molecular gas depletion times measured 
in many nearby spirals. 

ies at 15" resolution (corresponding to physical scales between 200 pc and lkpc) from 
Bigiel et al. (2010a) SFRs are estimated from GALEX far UV emission and S gas is es- 



timated from HI emission alone, assuming a negligible contribution from molecular gas 
on ^kpc scales in outer galaxy disks. For comparison, the filled contours show sampling 



data from the star forming disks of 7 spiral galaxies from Bigiel et al.(2008) (compare the 



right panel in Figure [5] above) . In the outer disks (open contours) S§fr scales with S gas , 
i.e., Shi, though the scatter in Ssfr for a given Hi column is large. This is an interesting 
difference compared to the inner parts of spiral galaxies, where the H 1 showed no clear 
correlation with Esfr (compare Section [2]). 

The Hi-FUV relation observed for outer disks suggests two things. First, that at large 
radii the availability of Hi may be a bottleneck for star formation. Even if stars form 
directly from H2, molecular clouds must be assembled from Hi and this will only be 
possible in regions with enough H 1 to assemble these clouds. Second, in the inner parts 
of galaxies many physical conditions important to the H 1-H2 conversion change while Shi 
remains approximately fixed, but in the outer parts Shi varies while other environmental 
conditions show comparatively little variation. As a result, Shi transitions from being a 
relatively unimportant driver for SF in the inner parts of galaxies to a key quantity at 
large radii. This is apparent from the right panel of Figure01 where we use extremely deep 
FUV data from GALEX and Hi data from THINGS to trace S F and Hi out to almost 
four optical radii in the nearby spiral M83 ( |Bigiel et al. (2010b) ). The Hi depletion time 
(H i-to-SFR ratio) inferred from the radial profiles in this plot is approximately constant 
at large radii: it is about a Hubble time, i.e., much longer than the ~ 2 Gyr molecular gas 
depletion time scale observed in the inner parts of galaxies. This suggests relatively fixed 
conditions for molecular cloud formation and that assembling H2, rather than forming 
stars out of H2, is the rate-limiting process for SF in outer disks. 
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4. The Composite Star Formation Law 

The combined distribution (inner and outer disks) in the left panel of Figure H] can 
be divided into different parts according to gas column density, each part describing the 
relation between gas and SFR in a particular regime in a typical spiral galaxy disk. For 
low gas columns (outer disks), the SFR scales with gas column, though with significant 
scatter. At smaller radii and higher gas columns — corresponding to much of the area 
inside the star-forming disk — the distribution becomes much steeper. In this regime, 
knowing S gas alone is not enough to predict Ssfr with any accuracy. Across this regime 
the H 2 -to-H I ratio is varying steadily as a function of other environmental quantities. 
At yet smaller radii and higher gas columns, the dominant phase of the ISM transitions 
from atomic to molecular and a strong correlation emerges between S gas and Ssfr- 



There is good observational evidence (e.g., Kennicutt(1998) Gao fc Solomon(2004) 
Greve et al.(2005)| Bouche et al.(2007) I that at higher gas columns, the relation steepens 



further, so that the SFR-per-H2 ratio is higher in starburst galaxies than in normal galaxy 
disks. This may drive the frequent observation of N ~ 1.5 in starburst galaxies. However, 
it must be emphasized that there is currently a lack of data probing normal disk galaxies 
and starbursts in a self-consistent way, so the details at the high end of this relation 
remain uncertain. 



5. Conclusions 

With vast improvements in the data available for nearby galaxies some consensus is 
beginning to emerge on how different parts of galaxies populate the SsFR-S gas parameter 
space. The role of environmental quantities other than gas surface density alone are 
beginning to become clear and different relations arc emerging for different types and 
parts of galaxies. When viewed in detail the composite relation may not be a simple power 
law, but it contains key information to constrain theories and to benchmark simulations. 

Challenges remain, too. The determination of star formation rates at low surface bright- 
ness is still difficult. The use of CO to trace H2 underpins almost all of this work but the 
CO-to-H2 conversion factor remains imprecisely calibrated as a function of environment. 
Finally, the fundamental units of star- formation, individual molecular clouds, remain 
largely observationally inaccessible beyond the Local Group — a situation that will not 
change until ALMA begins its full operations. 
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